Criteria to define a pair-ion plasma and 
the role of electrons in nonlinear dynamics 
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A criterion to define a pure pair-ion (PI) plasma is presented. It is suggested that the lighter 
elements (like H and He) are more suitable to produce PI plasmas. The observation of ion acoustic 
wave (IAW) in recent experiments with fullerene plasmas clearly indicates the presence of electrons 
in the system. A set of two coupled non-linear differential equations has been obtained for PI plasma 
dynamics. In moving frame, it can be reduced to a form similar to Hasegawa-Mima equation but it 
does not contain drift wave. Criteria, pair-ion plasma, ion acoustic wave 
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I. HISTORY OF THE PROBLEM 

Peculiar experimental observations [1] of pair-ion 
(PI) fullerene C$ plasma have invoked a great deal of 
interest in this topic. It has been reported that a pure 
PI fullerene plasma can support three kinds of electro- 
static waves propagating parallel to the external static 
magnetic field. These waves are the ion plasma wave 
(IPW), the ion acoustic wave (IAW), and the third one 
has been named as the intermediate frequency wave 
(IFW). In an earlier experiment [2], an alkali- metal- 
fullerene plasma (K + , e~ ,Cq ) was produced by intro- 
ducing fullerenes into the potassium plasma to realize 
a PI plasma. The IAW speed in Rcf. [1] has been 
defined as c s — where m + = m_ = m; and 

T + = T_ = Ti have been used. Here ji is the ratio 
of ion specific heats. The subscripts plus and minus 
denote the singly charged positive and negative ions, 
respectively. It is important to note that here c s k is 
the frequency of the ion thermal wave and is not the 
IAW. This thermal mode is similar to the sound wave 
in neutral fluids. 

The electrostatic waves were excited in fullerene 
plasma externally and following observations were 
noted[l]. First the IAW has frequency larger than the 
theoretically calculated frequency VTik where VTi — 
(£iS)3 is the ion thermal speed in our notation. Since 
electron density has been assumed to be zero (n eo = 0), 
therefore the ion acoustic speed, say c s , has not been 
defined as a function of electron temperature as c s = 

(S.)4. 

Second, it has been noticed that IFW has a feature 
that the group velocity is negative but the phase veloc- 
ity is positive i.e. the mode is like a backward wave. 
However, the IPW shows no special features in PI plas- 
mas. 

After these observations, some theoretical investi- 
gations have shown that the acoustic speed becomes 
larger in a pair-ion plasma if it is not pure and con- 
tains significant concentration of electrons [3,4]. The 



IAW and some other modes have been investigated in 
pair-ion-clcctron plasma and the linear IAW dispersion 
relation has been obtained using quasi- neutrality [3] . In 
Ref. [4] , the IAW in PIE plasma has been discussed us- 
ing kinetic approach and it has been pointed out that 
the quasi-neutrality is not a good approximation for 
PIE plasmas. As the number density of electrons de- 
creases, the electron Debye length Xdc — ( ^ nl e 1 ) 1 m ~ 
creases and the charge reparation effects become impor- 
tant. Furthermore, it has been shown that the Landau 
damping of IAW is reduced in the PI plasmas due to the 
presence of electrons and hence this mode can be ex- 
cited in the PIE plasmas easily in the limit 1 < Xjj e k 2 . 
The basic definition of plasma requires X^ e k 2 < 1, but 
the opposite limit k\ 2 De k 2 seem to be a possibel case 
in PIE plasmas. The experimental set up may be very 
reliable to produce PI plasmas, but the counter check is 
necessary to be sure whether the produced plasma can 
behave as a pure PI plasma or not. For this one needs 
to estimate the electron density n e o or the densities of 
positive n+o an d negative rt_o ions in the system. 

One cannot say with certainty that a plasma has 
exactly zero electron density. As the electron density 
is reduced in a system, the electron plasma frequency 
= ( 47r " e ° e may become smaller than the fre- 



ed 



quency of oscillations of positive ions viz, 



u pe 



« UJ 



(1) 



where, 



/47rn+o e ' : 



We need to find out a criterion for a PI plasma. It 
is the quantitative limit on the ratio of electron den- 
sity to positive ion density (^ £2 -) which can decide if 
the role of electrons is negligible in a plasma. Sev- 
eral authors [5,10] have theoretically studied various 
aspects of linear and nonlinear waves and instabilities 
in pure PI plasmas after the above mentioned experi- 
mental observations. In an investigation, the behavior 
of the so called backward mode (IFW) has been at- 
tributed to the nonlinear dynamics of PI plasma [5]. 
Recently two more experimental research papers have 



appeared in the literature on PI plasmas [11,12]. In 
one of these works [11], an effort has been made to 
produce PI plasma with positive and negative Hydro- 
gen ions (H + + H~). On the other hand, already a 
recent theoretical work [13], it has been pointed out 
that a plasma of lighter ions (of low - Z materials) can 
behave as a pure PI plasma with relative larger elec- 
tron density n e o compared to plasmas of heavier ions 
like fullerenes. Therefore, it has been suggested that 
it is more suitable to try to produce PI plasmas us- 
ing low-Z materials if other physical conditions like the 
electron attachment cross-section (to produce negative 
ions), ionization and recombination rates can be con- 
trolled. Some linear and nonlinear waves of PI and PIE 
plasmas are also discussed in this work. 



II. IAW IN PIE PLASMAS 



where N = ( " + "+ "~° )- This is the same as Eq. [5] 
of Ref.[4]. Note that c s k < uj s when n e0 < n +a because 
1 < N . The IAW can be excited easily in PIE plasma 
because the Landau damping rate decreases in the limit 



III. CRITERIA FOR PAIR-ION PLASMA 

To decide whether the produced plasma can be 
called a pure PI plasma or not, one needs to estimate 
the ratio ^— ; . If the electrons pressure is not unusually 
high due to some external heating mechanism, then 
the situation uj 2 e « u 2 i+ implies that their role in 
plasma dynamics can be neglected. Such a system 
can behave as a pure PI plasma. Let us consider the 
plasma wave dispersion relation using kinetic approach 
written as, 



Here we derive the linear dispersion relation of ion 
acoustic wave (IAW) in pair-ion-electron (PIE) plasma 
for the simplest case, T + = T_ = T» «T e using fluid 
equations. The wave is assumed to propagate along 
the external magnetic field to compare the theoretical 
result with the experimental observation [1]. 

Let the constant external magnetic field be along z- 
axis, i.e. B = and the wave be propagating along 
the field lines with wave vector k"*= kz. The equations 
of motion for singly charged positive and negative ions 
become, respectively, 



and 



d t v +z ~ d z ip 

mi 



d t v- z ~ H d z ip 

rrii 



(2) 



(3) 



E 



{1 + L^ZjW(Zj)} = (8) 



where Af^- 



iminj e 2 ' i ~~ sp2kv T j 



V Tj = and 

W(Zj) is the plasma dispersion function for the jth 
species (J = ±,e). In equilibrium, the quasincutrality 
demands n 0e + n_ = n +0 where both the positive and 
negative ions are assumed to be singly charged. This 
is the case of Ref. [1] as well. For the time being, 



let us assume Tj < Te and vx e « 



such that 



1 < < \Z e \ < < \Zi \ which allows us to use an asymptotic 
expansion of W(Zj) to study the system analytically. 
The case of IAW with vt% « j « «Te has been 
discussed in Ref. [4] in detail. 

Assuming m+ = m_ = m,, Eq. [8] can be expressed 
as, 



The continuity equations yield, 

d t (n + - rc_) + n +0 d z v +z - n^d z v^ z = (4) 

Assuming electrons to follow the Boltzmann density 
distribution in electrostatic field E — —Ap as, 



W --P>pi+ ^ k VTeU; e +k 2 vi i -u; i _+k 2 vi i+ +uZ i+ } 



2 J Z e 



Zi. 



De 



_ P -Zi + Zl+ P -z i+ 

K A Di+ 



n e ~ n e0 e T » ~ n e0 (l + — ) 



(5) 



(9) 



the Poisson equation gives, 



(n+ -n_) 



(6) 



Using the above set of equations, we can obtain the 
linear dispersion relation as, 



where P° = (1+2=2 + ™e<L so.) . For u = w r -i7, the real 

n v ro+o ™+o m e ' 1 

and imaginary parts of Eq. [9] become, respectively, 



2 ^ D 2 ^r2 i 2 I / 10 pi- i 2 \ 2 \ 2 

UJ r — f n U) pi ± + — q [V Ti ±K + (—2 K ^Di-i^pi- 

r„ UJ pi+ 



u 2 s = iVo 



l + \ 2 De k 2 



(7) 



+ (-^k 2 \ 2 De )u : 



pel 



(10) 



and 



plasma: 



Z, 



i+ 



exp{~Z? + } + 



Zj- 



+ 



tc 2 \ 2 



sxp{-Z 2 }]w r 



■exp{-Zf_} 



(11) 



The important point to note is that in the limit to pe < 
uj p i + and VTek < w the electron plasma wave turns 
into the ion plasma wave. Furthermore in the limit 
(n 2 ^) ~~ ^' we nee( l not to use the perturbed electron 
density n e i in the Poisson equation. 

The comparison between the two terms of Eq. [10], 

, can be very important to 



J Te 



ui 2 e and Art4j_| 



decide if the plasma can be called a pure (PI) system. 
Note that if Tj = Ti + = and = = rrij, 
Ti w pi+ provided that 



then k 2 v^ e LOp e « k 2 v%,J 1 



n 0e <£ . T i ( m e 
n +0 " T e \m l 



(12) 



In this case the electron contribution to the plasma 
dynamics can be neglected only for a very small value 

of the ratio . Then the plasma can be called a pure 

««+ 

(PI) plasma. 

In the case of fullerene plasma, = 720m p (where 
m p is mass of the proton) and ^ <~ y^gg- Therefore 

we have ^ = 7.56 x 1(T 7 . If T; < T e is assumed, then 
the fullerene plasma discussed in Ref. [1] can be called 
a pure pair-ion plasma only if the following limit holds: 



^ « (7.56 x 10" 7 ) 2 ^ 



(13) 



c s k « akv T e « u pi+ 



(16) 



If we choose m e /rrii « a « 1 the condition (16) is 
satisfied and it is then in agreement with the fact that 
the ion acoustic wave should not appear in the pure (PI) 
plasma. However, a smaller value of a is preferable. If 
n eo is so small that a is almost zero, then we will have, 



av 2 Te k 2 «c 2 s k 2 <«lu 



(17) 



In this case the IAW remains almost non-existent. 
That is n eo is too small and electron pressure does not 
contribute to plasma dynamics. The only normal mode 
of the system with k||B is the ion plasma wave which 
may have a negligible contribution from the small num- 
ber of hot electrons. For the case of Helium (He) plasma 
m e /m p ~ 10~ 4 and if T, < T e is assumed, then it will 
become almost a pure (PI) plasma if noe/n-i « 1CT 5 
holds. Therefore we conclude that it is more suitable 
to try to produce (PI) plasma of lighter atoms (or 
molecules) if other physical conditions like the ioniza- 
tion/recombination rates, and the electron attachment 
cross section can be controlled. Therefore it is sug- 
gested that the Hydrogen and Helium systems can be 
very useful to achieve a pure PI plasma. 



IV. VORTICES IN PIE PLASMAS 

If Ti ^ 0, then the perpendicular drift velocities for 
ions can be written as, 

Vbi X Z 1 



V,_L 



B 



-Ei x z - 



- 7T-(Vi + Vj.V)vj x z 



The plasma density in Ref. [1] is m ~ 10 7 cto~ 3 . It 
means that this system can become a pure (PI) plasma 
only if there is no electron in the system which is very 
unlikely physically. 

Fortunately, we have a better condition than (12) 
to call the plasma a pure (PI) plasma. Our main re- 
quirement is the limit oj pe « uj±, which replaces the 
relation (12) by a new limit as follows: 



n 0e 

n 



« a- 



rm 



(14) 



V£ + V D j + V. 



P3 



(18) 



Here j = ± and fL = p or electrons, we have 

c V„i x z „ „ 

Vel = — E ± x z - = V B + V De 19 

Bq il e m e n e 

where \d t \ « O e = has been used. The continuity 
equations of the ions yield, 



where a must satisfy the condition a << 1. Cor- 
respondingly the condition on thermal correction term 
is, 



c c 
d t (n + - rc_) + — Vn e0 .(Z x V_l ) - -5-7^(^+0 + "-0) 



k 2 v 2 



P e <• , ,2 

Te —2— ^ W pi+ 



(15) 



Since 



a therefore (15) suggests the following 



maximum limit on the value of a to call a plasma as PI 



x (d t + v B .V)V 2 = n-o<9z v z - - n +Q di v z+ (20) 
and the parallel equation of motion becomes, 
(dt + Ve-V) = n-od z v z - - n +0 diV z+ 



= — (n+o + n- )d z (j) 
rrii 



(21) 



V. NONLINEAR EQUATIONS FOR PI 
PLASMA DYNAMICS 



Assuming Boltzmann density distribution for electrons 
n e <~ n e t)e T ' and using the Poisson equation. 

V.{E) = 47re(n+ - n_ - n e ) (22) 

the nonlinear Eqs. [19] and [20] can be written, respec- 
tively, as: 

dt{~\ 2 Dl N 2 <S>} + D e K ne .(z x Vi$) - N p 2 s 



and 



where V 
D 



x(d t + D e zx V_L$.V)Vi$ = 9 Z F (23) 



(d t + D e z x V ± $.V)F = c 2 s N d z $ (24) 



_2 

and c s = 



__ (ra +0 +w-n) _ e0 



= m anci c * = fe) 3 - 

Equation [23] is the Hasegawa-Mima (HM) equation 
for (PIE) plasma if the RHS is ignored (for v iz -> 0) 
and X 2 De k 2 « 1 is assumed. 

These equations give the coupled linear dispersion 
relation of drift wave and IAW in (PIE) plasmas as, 



G oj 2 



N c 2 s kl = 



(25) 



where G = (1 + X 2 De k 2 + N p 2 k 2 j_) and uj* e = vg.k, 

If N c 2 k\ «w* holds, then we obtain only the drift 
wave dispersion relation as, 



(26) 



(1 + X% e k 2 +N p 2 kl) 



In (PIE) plasma the quasi-neutrality can break down 
for IAW in the limit 1 << \ 2 De k 2 because n e o can be 
very small. On the other hand the inequality Xo e < Ps 
always holds. Since p 2 k 2 < 1 in the fluid model, 
therefore X? De k 2 should not be much larger than 1. 
This means in magnetized plasmas, the IAW cannot 
have wavelengths shorter than Xoe within fluid theory 
framework because of the limit, 



*Dek' Z < P\k\ < 1 



(27) 



It is important to note that as n e Q decreases, -/Vg in- 
creases to have \ 2 De k 2 « N plk"j_ and hence Eqs. [23] 
and [24] give the pair plasma convective cell (PPCC) 
mode, 



As n en decreases, the IAW converts into the PPCC 
mode. In between these two limits, the electron drift 
wave couples with IAW and PPCC. 



The set of nonlinear Eqs. [23] and [24] can be trans- 
formed into HM equation in a moving frame which ad- 
mits monopolar and dipolar vortex solutions [13]. The 
nonlinear dynamics of (PI) plasmas are described by 
Eqs. [23] and [24] in the limit 1 << 7V and they re- 
duce, respectively, to the following equations: 

(dt + DiZX Vi$.V)Vi$ = --^^y (28) 



and 



(d t + Diz x V ± $.V)V = 2v 2 Ti d z <§> 



(29) 



where $ = §£, A = ^ Pi = ^ and V = (v z --v z+ ) 

In the stationary (77, x) frame, the coupled Eqs. [28] 
and [29] can be written as, 



Cid„Vi$ + Vod n $ + {Vi$, $} = 



(30) 



where Ci = ^, V = ^L 2p 2 D u L - \ C - 

and Co is an arbitrary constant. The important point 
to note is that the form of Eq. (28) is similar to 
Hasegawa-Mima equation but the physics of the equa- 
tion is completely different. The set of nonlinear Eqs. 
[28, 29] is valid as well for electron-position plasmas in 
the classical limit. But these equations do not contain 
the drift wave and the ion acoustic mode. They de- 
scribe the nonlinear dynamics of (PI) plasmas in the 
quasi-neutrality approximation. 



VI. PI PLASMA FOR FUSION 



The drift waves are the fundamental source of insta- 
bilities in tokamak fusion plasmas. Several kinds of re- 
active and dissipative drift instabilities appear in Toka- 
mak plasmas. The presence of electrons in laser-fusion 
is also problematic for laser absorption, heat conduc- 
tions and uniform compression. The Pi-plasmas can be 
very suitable fuel for fusion, in principle because drift 
waves cannot exits in such plasmas. However, it does 
not seem easy, at least at present times, first to produce 
PI plasma of Hydrogen or Helium at high densities and 
high temperatures for fusion. Second the confinement 
for a longer time can also be a problem because of re- 
combination and production of electrons and neutrals 
as a result of collisions between positive and negative 
ions as well as with neutrals. 



VII. SUMMARY 

A criteria to define a pure pair-ion (PI) plasma has 
been discussed. The condition (14) must be satisfied 
along with (16) to call a plasma as a pure PI plasma. 
But a very small value of a such that a « — is 
preferable. It is suggested that the lighter elements 
are preferable to produce pure PI plasma if electron 
attachment cross-section and recombination rate can be 
controlled for desirable results. In the plasma of lighter 
elements, the condition (14) can be satisfied even for 
relatively larger values of a. 

It has also been stressed that the observation of ion 



acoustic wave (IAW) in the experiment, itself is an in- 
dication of the existence of significant concentration 
of electrons in the produced pair-ion fullerene plasma 
and hence it cannot behave as a pure PI plasma. Fur- 
thermore, the frequency of IAW increases in a pair-ion 
plasma in the presence of electrons because we have 
1 << N . 

A set of two coupled nonlinear differential equations 
has also been obtained for the pure PI plasma. In a 
moving frame, these equations can be reduced to a sin- 
gle equation similar to Hasegawa-Mima equation but it 
does not contain drift wave. 
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